1 



Long-Range Spectral Statistics of Classically Integrable Systems 

Investigation along the Line of the Berry-Robnik Approach 

Hironori MAKING ■^'*^ and Shuichi Tasaki^ 

^Department of Human and Information Science, Tokai University, Hiratsuka 

259-1292, lapan 

'^Department of Applied Physics, Waseda University, Tokyo 169-8555, lapan 

(Received April 25, 2005) 

Extending the argument of Ref^ to the long-range spectral statistics of classically inte- 
grable quantum systems, we examine the level number variance, spectral rigidity and two- 
level cluster function. These observables are obtained bvapplying the approach of Berry 
and RobnilP and the mathematical framework of Pandey-* to systems with infinitely many 
components, and they are parameterized by a single function c, where c = corresponds to 
Poisson statistics, and c 7^ indicates deviations from Poisson statistics. This implies that 
even when the spectral components are statistically independent, non-Poissonian spectral 
statistics are possible. 

§1. Introduction 

For bounded quantum systems in the semiclassical limit, the statistical prop- 
erties of eigenenergy levels have been studied intensively. Universal behavior has 
been found in the statistics of unfolded energy levels which are the sequences of 
numbers uniquely determined by the mean level density in a small energy interval 
obtained using the Thomas- Fermi rule. It is widely known that for quantum systems 
whose classical counterparts are strongly chaotic, the energy level statistics are well 
characterized by random matrix theory (RMT), in accordance with the conjecture 
of Bohigas et al., ' which gives a Wigner-type level-spacing distribution. For sys- 
tems whose classical counterparts are integrable, the energy level statistics are well 
characterized by Poisson statistics in accordance with the conjecture of Berry and 
Tabor,!^ which gives a Poisson-type level-spacing distribution. 

Universal behaviors are also exhibited by other quantities, e.g., the level number 
variance i7^(L), which is the average variance of the number of levels in an interval 
containing an average of L levels.!^ For an unfolded level sequence, this interval is 
equivalent to an interval of length L, and ^'^{L) is thus defined as 

S\L) = {{MiE + L) - MiE) - Lf), (M) 

where J\f{E) is the number of eigenenergy levels below E, and the brackets (• • • ) 
represent the average over the value of E. The number variance depends primar- 
ily on the long-range correlation between the energy eigenvalues, in contrast to 
the level spacing distribution, which depends primarily on the short-range corre- 
lation. In RMT, for which /? = 1,2, and 4 correspond, respectively, to GOE, GUE, 
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and GSE level statistics, the number variance increases logarithmically with L 
[explicitly, S'^{L) ~ (2//37r^) log(27rL)] , while in Poisson statistics, the number vari- 
ance is equal to the number itself [i.e., S'^{L) = L] . The latter behavior follows from 
the fact thatZ'2(L) can be expressed in terms of the two-level cluster function^' 12 (r) 
as 

E'^{L)=L-2[ {L - r)Y2{r)dr, (1-2) 
Jo 

where Y2{r) is ordinarily defined in terms of the two-level correlation function.!^ 
In the absence of spectral correlations, we have Y2{r) = for all r, which implies 
E^{L) = L. 

The level statistics for integrable quantum systems have been theoretically stud- 
ied in a number of works,^''^*^''^"*^'' i^''^'^'' i^'' I5i,i i,i6i, 17 1, 18, i ^^-^^ have been the subject 
of many numerical investigations. In RefsllOj) and I17j) . eigenvalue problems for in- 
tegrable quantum systems are reformulated as lattice point problems. There, it is 
proved, under explicit diophantine conditions, that the two-level correlation function 
exhibits Poisson statistics. 

The appearance of Poisson statistics is now widely acknowledged to be a univer- 
sal phenomenon in generic integrable quantum systems. However, the mechanism 
responsible for this statistics is still not well understood, and deviations from Poisson 
statistics have also been found.El'l^'E^naJiniCsi.iei.ui 

As suggested, e.g., by Hannay (see the discussion given in Ref,2)), one possible 
explanation for the appearance of Poisson statistics is the following: For an integrable 
system of / degrees of freedom, almost every orbit is generically confined in its 
inherent torus, and the entire region in the phase space is densely covered by invariant 
tori, as suggested by the Liouville- Arnold theorem.-^ ' In other words, the phase space 
of an integrable system consists of infinitely many tori that have infinitesimal volumes 
with respect to the Liouville measure. Thus, the energy level sequence for the whole 
system consists of the superposition of the subsequences that are contributed from 
those regions. Therefore, if the mean level spacing of each independent subset is 
large, one would expect Poisson statistics as a result of the law of small numbers.!^ 
This explanation given by Hannay is supported by the theory proposed by Berry 
and Robnik.EJ 

The Berry- Robnik theory relates the formation of the eigenenergy levels to the 
phase-space geometry by assuming that the sequence of the entire energy spectrum 
is given by the superposition of statistically independent subspectra, which are con- 
tributed respectively from eigenstates localized in the invariant regions in phase 
space. The existence of such independent subspectra is a consequence of the con- 
densation of energy eigenfunctions on disjoint regions in the classical phase space 
and the resulting lack of overlap among these eigenfunctions. Therefore, indepen- 
dent subspectra can be expected only in the semiclassical limit, in which the Planck 
constant tends to zero (/i — > 0). This mechanism was first proposed by Li and Rob- 
nik^^' on the basis of an implicit state described by Berry,^^''— ^^d it is sometimes 
referred to as the principle of uniform semiclassical condensation of eigenstates^^ 
This principle states that the Wigner function of a semiclassical eigenstate is con- 
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nected on a region in phase space explored by a typical trajectory of the classical 
dynamical system. In integrable systems, the phase space is folded into invariant 
tori, and the Wigner functions of the eigenstates tend to delta functions on these tori 
in the semiclassical limitP^ On the other hand, in a strongly chaotic system, almost 
all trajectories cover the energy shell uniformly, and hence the Wigner functions of 
eigenstates are expected to become delta functions on the energy shell, as suggested 
by the quantum ergodicity theoremP^''^ Because of the suppression of tunneling, 
each quantum eigenstate is folded into independent subsets in the semiclassical limit, 
and is expected to form independent spectral components. Indeed, the formation of 
such independent components has been confirmed numerically in a deep semiclassical 
regime.ESJ^ 

Employing the above point of view, Seligman and Verbaarschot investigated 
a semiclassical formula for the number variance S^iL) in the case that the en- 
ergy spectrum consists of the superposition of subspectra contributed by the phase 
space components.!^ They introduced the mathematical framework formulated by 
PandejJ^ and observed that the long-range spectral statistics of systems consisting of 
N independent components described by GOE statistics of equal strength approach 
Poisson statistics in the limit N — > +00. 

The level statistics of systems with infinitely many spectral components have 
also been studied in Refs. [71) ,121) and i3l|) . For example, in Ref. [T|), the level 
spacing distribution in the limit A'^ — > +00 is investigated, and it is classified into 
three cases: (i) Poisson distribution, (ii) Poisson distribution for large level spacings, 
but possibly not for small spacings, and (iii) sub-Poisson distribution. It is observed 
that cases (ii) and (iii) are possible when the spectral components exhibit strong 
accumulation and have a singular statistical property. 

These works seem to imply that the existence of infinitely many independent 
components is an essential ingredient for the appearance of Poisson statistics and 
that the semiclassical limit is one situation in which an effectively infinite number 
of components arises. However, in some classically integrable systems in which one 
might expect there to be infinitely many components, deviation from Poisson statis- 
tics has been observed .^''^'^'^''^'^'^'^ Therefore, it is interesting to explore 
the consequences of the existence of infinitely many independent components. 

Here we briefly review the formula proposed by Pandey.!^ Consider a system 
whose classical phase space is decomposed into A^-disjoint regions which provide 
independent spectral components. When the entire level sequence is a product of a 
statistically independent superposition of subsequences, the level number variance 
for the entire sequence, ©Z"^, can be decomposed into those of the subsequences, Sf, 
as 

N 

U\L;N) = J2^f{PiL), (1-3) 

1=1 

where pi{i = 1,2,3, ■■ ■ , N) denotes the statistical weights of the individual spectral 
components. This quantity is defined by the mean level spacing Si of the unfolded 
spectral components as pi = 1/ Si and satisfies X^^i Pi = 1. In terms of the cluster 
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function Y2^i{r) of a subsequence, Uf{piL) is given by 

Ef{p,L) =piL-2 r^(/9.L - r)Y24r)dr, (1-4) 
Jo 

where satisfies the relation Y2^i{r)dr = 1. 

In combination with the number variance, the spectral rigidity ^^{L) introduced 
by Dyson and Mehta has played a major role in the study of long-range spectral 
statistics!^ It is defined as the mean-square deviation of N{E) from the best-fitting 
straight line; 

I / i-Es+L \ 

Z\3(L) = - ( min / dEiNiE) - AE - Bf) . (1-5) 

L \ A,B Je^ I 

The angular brackets (• • • ) represent the average over the value of Eg. Because the 
number variance is connected to the spectral rigidity aP 

Az{V) = dr{L' - 2L\ + r^)E^ir), (1-6) 

A^lL) is also additive and the formula (|l-3p can be rewritten as 

TV 

As{L;N)=J2As,i{piL), (1-7) 
1=1 

where the rigidity of the subspectrum is described by the number variance of the 
subspectrum as 

In the argument of Berry and Robnik,^^ the statistical weights of the individual com- 
ponents in the formulae l\l-'6\i and (|1-7|) are equivalent to the phase volumes (Liouville 
measures) of the corresponding invariant regions. This relation is satisfactory in a 
given small energy interval if the Thomas-Fermi rule for the individual phase space 
regions still holdsP Therefore, the formulae and (|l-7() relate the level statistics 
in the semiclassical limit to the phase-space geometry. 

In most general cases, the cluster function might be singular. In such a case, it 
is convenient to use its cumulative distribution function cf. 

Ci{L-pi)=2 r Y2,iir)dr. (1-9) 



In addition to Eqs. and (|1-4|) . we assume the following two conditions for the 

statistical weights: 

• Assumption (i): The statistical weights of individual components uniformly 
vanish in the limit of infinitely many regions: 

max/jj^O as ^ +oo. (I'lO) 

i 
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• Assumption (ii): The weighted mean of the cumulative cluster function, 



N 



c{r;N) = '^piCi{r;pi) 



(Ml) 



converges in the +00 limit to c(r): 



lim c(r; A^) = c(r). 



(M2) 



The limit is uniform on each closed interval in the range < r < L. 
Under Assumptions (i) and (ii), formula (|1-3|) leads to the following formulae in 
the limit N ^ +00: 



When the mean level spacings of individual components are sufficiently sparse {i.e., Si - 
+00), we expect that c(r) = 0, and the limiting observables (|1-14|) and (|1-14|) of 
the whole energy sequence reduce to those of Poisson statistics, ^|=o ~ ^^'^ 
As,c=o = L/15. In general, we expect c(r) 7^ 0, which corresponds to a certain 
clustering (accumulation) of the levels of individual components. 

In the next section, the above statement is proved, and ^'^{L; N) and ^^{L; N) 
in the limit A^ — > +00 are also investigated. In section |3J the cluster function Y2 of 
the whole energy levels is also analyzed in the A^ — > +00 limit. We give one possible 
example of the deviation from Poisson statistics in section^ and in the concluding 
section, we discuss some relationships between our results and those of related works. 



In this section, starting from Eqs. (|l-3() and ()l-4() and Assumptions (i) and (ii) 
introduced in the previous section, we show that, in the limit of infinitely many 
components (A^ —>■ +00), U'^{L;N) and Z\3(L;A^) converge to the formulae (|l-14l) 
and Hl-14() . respectively, with the cumulative cluster function c. The convergence is 
shown as follows. 

First we rewrite U'^{L;N) in terms of the cluster function l2,i of independent 
components: 




(1-13) 



(1-14) 



§2. Number variance and spectral rigidity 




(2-1) 



Now, note that using Eq. (|l-9j) and integration by parts, we have 




(2-2) 
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Since the eonvergence X^^i piCi{r; pi) c(r) as N —>■ +00 is uniform on each interval 
r e [0, L] by Assumption (ii), 17^ (L; N) has the following limit: 



U\L; N) = L- [ V piCi{r; pi)dr (2-3) 

•^0 ^=l 

— > L- j c{r)dr as N ^ +00. (2-4) 
Jo 



'0 

Therefore, we obtain the desired results: 



hm U\L- N) = El{L) = L- T c{r)dr, (2-5) 

A^^L) = ^ driL^ - 2LV + r^) lim E\r; N) 

= T^ + TaI {L^ - + r^) [ dxc{x). (2-6) 
15 -L* Jq Jo 

We have Poisson statistics, i.e., = L and /^3,c = L/15, if the energy levels are 
uncorrelatcd, and thus c(r) = for all r. Such a case is expected when the individual 
cumulative function Ci{r;pi) is bounded by a finite positive function g{r) as 

\ci{r-pi)\<p1g{r), (2-7) 

for all i and > 0. Indeed, we have 



N N 



|c(r;iV)| < |ci(r;pi)| <5(r)^PiP,'' <5(r)(maxpif ^0 = c(r), as N 



1=1 



More specifically, for example, we have c = if the absolute values of individual 
cluster functions are uniformly bounded by a positive constant D: \Y2^i{r)\ < D 
{I <i < N). Indeed, the following holds: 



\<L-N)\ 



•^^Pi r'^Y2,i{r)dr 
i=i Jo 

N 

< 2DL pf < 2DL max: pi^^pi — >0 = c{L) , as AT ^ +o(^2-9) 



N N 



i=l 



This includes the case in which the overall energy spectrum consists of the super- 
position of N chaotic components characterized by GUE level statistics, Y2^i{r) = 

^ smjTrr) j . ThuS, WC haVC 
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N 

1=1 

N 

< max /f3i Pi ^ as ^ +00, (2-10) 

i=l 

where we have used the relation sino" < a for a >0 and Assumption (i). 

In general, we expect that c(L) / 0, which corresponds to a strong clustering 
of the levels of individual components, leading to a singular cluster function Y2^i{L) 
and a nonsmooth cumulative function Ci{L;pi). Because of the singularity of l2,j 
in Eq. ()l-9() . Cj might converge to a nonzero constant even in the limit N ^ +00. 
For example, if all Cj(r) simultaneously satisfy Cj(r) > C > or Cj(r) < —C < at 
some two-level distance r with a positive constant C, there is deviation from Poisson 
statistics, and we have |c(r)| = limjv^+oo X^i^i PilciC?^)! > C > 0. Such a strong 
clustering would be expected when there is a symmetry in the system. One example 
of a system that might possess such a symmetry is a rectangular billiard system. 
In section [IJ we show that such systems exhibit non-Poisson level statistics with 
c(L) ^ 0. 



§3. Two-level cluster function 



In this section, we discuss the two-level cluster function of the entire level se- 
quence, 

Y2{L-N) = -^U\L;N), (3-1) 

in the limit of infinitely many independent components, +00. 

Substituting Eqs. (|l-3p and (|1-4|) into the above equation, Pandey's formula 
Hl-3|) can be rewritten for the cluster function as 

j2 ^ AT 

Y2{L;N) = -—^uf{p,L) = 2 ^ p2^2,.(p2^). (3-2) 
1=1 1=1 

The second equation follows from the commutative relation between the sum and 
the differentiations. However, in the limit +00, this relation does not generally 
hold, due to singularities of the sum. Therefore, in the present work, we investigate 
the limit of the two-level cluster function only for the case of the weak convergence 
to Y2^ciL), i.e., the case in which 

Y2AL) = -^^c- (3-3) 



According to Helly's theorem,'23'E3 the weak convergence limit is defined by 

fL rL 

W->+oo 



lim / Y2{r;N)dr= [ Y2c{r)dr. (3-4) 
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This equation is equivalent to limjsf^^oo -^^"^{L; N) = ^i7|(L), and it can be 
proved by using Eqs. (|2-1() and (|2-2() and Assumption (ii) as follows: 



dL ' ' ' dL 

i=l 
N 

= I PiCi{L] Pi) (3-5) 

i=l 

— >l-c{L) = -^Sl{L) as N^+oo. (3-6) 

When the limiting function c(L) is differentiable, the asymptotic cluster function is 
given as follows: 

y2,-c{L) = ^c{L) (3-7) 
§4. Rectangular billiard systems 

As an example of deviation from Poisson statistics, we examine a rectangular 
billiard system whose energy levels are given by 

en,m = + a , (4-1) 

where n and m are positive integers, and a is the aspect ratio of the rectangle. The 
unfolding transformation of the eigenenergy levels {en,m} — > {^n,m} is carried out 
by using the leading Weyl term of the cumulative mean number AA(e) of the energy 
levels, 

J\f{en,m) = -^^en,m- (4-2) 

'a 



It is well known that the two-point correlation function exhibits Poisson statistics 
when a is an irrational number and satisfies the diophantine condition. This fact 
was first pointed out by Sarnacl/^ and was later proved by Eskin et al. in Ref . ITTjl . 
While the two-point correlation function is unbounded and does not exhibit Poisson 
statistics when a is a rational number, expressed as a = p/q, where p and q are 
coprime positive integers.^'^ Thus it is interesting to study the spectral statistics 
for rational values of a. 

Fig. ^ displays the numerical results of the number variance U'^ for rational and 
irrational values of a. The numerical computation was carried out using a double 
precision real number operation. In the case that a is an irrational number, e.g., 
the golden ratio, a = (1 + -^/5)/2, the number variance exhibits Poisson statistics: 
Z'^(L) = L [see Fig. ^a)]. By contrast, in the case that a is a rational number, e.g., 
a = 1, a strong clustering of levels appears, due to the corresponding geometrical 
symmetry of the system, and S'^{L) clearly deviates from Poisson statistics [see Fig. 

mb)]. 

Fig. [21 displays the numerical results of c(L) for the values of a corresponding 
to Figs, ^a) and (b), respectively. For each value of a, c(L) is calculated from the 
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overall spectrum as 



c(L) = 1 



(4-3) 



From this, it is clearly observed that c(L) = when i^^(L) exhibits Poisson statistics, 
while c(L) 7^ when S'^{L) deviates from Poisson statistics. 

In order to characterize the deviation from Poisson statistics, i.e., c(L) ^ 
behavior, in terms of the singular statistical properties of the spectral components, we 
identify each energy level with a component as follows. Consider a set of eigenstates 
in the (re, m^fa) plane, as shown in Fig. EJ where the energy interval [e, e + Z\e] is 
represented by the striped region. This interval is then divided into IN subintervals, 
each corresponds to an angle of AQ = ir/iN. The subintervals are numbered from 
1 to beginning in the middle of the interval and moving both clockwise and 
counterclockwise. Here, there are two subintervals corresponding to each number. 
Then, the energy level e^.m £ [^j e + belonging to each of the subintervals is 
identified with the i-th spectral component as 



AN . . 

arctan I — Wa I — N 

vr V re 



+ 1, 



(4-4) 



where [x] represents the maximum integer that does not exceed x. With this identi- 
fication rule, the relative weight of each spectral component, pi (i = 1, 2, 3, ■ • • , A^), 
is given by 

Pr = ^, (4-5) 

where the number of components, N, is given in terms of the constant values 7 = 
Ae/e = Ae/e and rej (the number of levels identified with i-th component) as 



N 



7 



(4-6) 



The limit of infinitely many components, N — > +00, corresponds to the high energy 
limit e — > +00, which is equivalent to the semiclassical limit. In this limit, we have 
A6 0, and pi satisfies Assumption (i). 

Figs. IHa) and (b) display the numerical results of c(L) = YliLi PiCi{L) for the 
values of a corresponding to Figs, and (b), respectively. In each case, Ci{L) is 
calculated from the number variance, Ef, of the spectral component as follows: 



c»(L) = 1 



1 d 

Pi dL 



UfipiL). 



(4-7) 



From this, it is clearly seen that c{L) ^ when a = 1 and c{L) = when a = 
(1 + Y^)/2. However, as shown in Fig. E^b), c{L) does not reproduce c(L) well, i.e., 
the deviation from Poisson statistics is not as large as that found for the value of 
c(L) obtained from the overall spectrum. This disagreement is introduced artificially 
by the identification method employed in this section. In this identification rule, a 
pair of energy levels accumulated strongly is decomposed into different components. 
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For this reason, the singularity characterizing the statistical properties of the whole 
spectrum is suppressed, and the contribution to this singularity from the accumulated 
pairs is limited only to the contribution from the pair of levels which are identified 
with the same component. 



In this paper, starting from the formula ()l-3() . we investigated the energy level 
statistics of bounded quantum system with infinitely many components and stud- 
ied its deviations from Poisson statistics. In the semiclassical limit, reflecting in- 
finitely fine phase space structures of the classical dynamical system, individual 
energy eigenfunctions are expected to be localized in the phase space and to con- 
tribute independently to the level statistics. With this in mind, we considered the 
situation in which the system consists of infinitely many components, each of which 
gives an infinitesimal contribution. Then by applying the arguments of Pande}^ 
and Berry and Robnik,^' the level number variance, the spectral rigidity and the 
two-level cluster function were obtained within the limit of infinitely many spectral 
components. These limiting observables are described by a single function, c(r), of 
the energy interval r, which is characterized by c(r) = for the Poisson statistics and 
by c(r) 7^ when there are deviations from Poisson statistics. Thus, even when the 
energy levels of individual components are statistically independent, non-Poissonian 
level statistics are possible. 

Note that a singular two-level cluster function can be taken into account by 
a non-smooth cumulative function q. Such a singularity is expected when there 
is strong clustering (accumulation) of the energy levels of individual components. 
One possible example containing such strong clustering is the rectangular billiard 
system. As shown in the section ^ this system exhibits deviations from Poisson 
statistics when the aspect ratio of the rectangle is a rational number. To obtain a 
suitable method of dividing each energy level into the components, it is necessary 
to carry out further studies of this system. Other examples have been studied by 
Shnirelman,12i^ Chirikov and Shepelyansky,'^ and Frahm and Shepelyanskj-^ for 
certain types of systems containing quasi-degeneracy of levels resulting from time- 
reversal symmetry. As is well known, the existence of quasi-degeneracy leads to 
non-Poissonian level statistics. 

It is also interesting to extend the Seligman-Verbaarschot formulcP^ for the level 
statistics of nearly integrable quantum systems with two degrees of freedom. The 
classical phase space of a nearly integrable quantum system consists of regular and 
chaotic regions. For this kind of system, Seligman and Verbaarschot proposed the 
number variance of the entire energy spectrum consisting of a subspectrum (i = 0) 
exhibiting Poisson statistics as well as subspectra (i = 1,2, 3, ••• ,N) respectively 
exhibiting RMT: 



§5. Discussion and conclusion 





i=l 
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where ^Poisson(-^) = ^^'^ ^«\rmt number variance of the ith chaotic component 

obtained from the RMT. Because the Liouville measure of a chaotic region is larger 
than zero, i.e., pi > 0{i > 1), Assumption (i) is invahd for this system. However, 
the regular regions consist of infinitely many subsets, and the present approach is 
partially applicable to the spectral components corresponding to regular regions. 
Along the line of the approach employed in section [21 the original proposal for the 
level number variance and the spectral rigidity by Seligman and Verbaarschot should 
be replaced by 

/ fPoL \ Ji^ 

SliL) = i^poL - c{r)drj + ^ IJl^^^{p^L). (5-2) 

When c(r) = for all r, the new formula H5-2|) is reduced to that of the original pro- 
posal by Seligman and Verbaarschot fFormula l5-l() . In the case that c(r) ^ 0, on the 
other hand, the formula (|5-2() allows deviations from the Seligman- Verbaarschot for- 
mula. Therefore, this result might represent an extension of the Seligman- Verbaarschot 
formalism proposed for nearly integrable quantum systems. 

Acknowledgements 

The authors would like to thank Professor N. Minami, Dr. Y. Sakamoto, Pro- 
fessor M. Robnik, and Professor A. Shudo for penetrating comments and discussion 
which helped the authors to deepen and clarify their own arguments. The authors 
also thank the Yukawa Institute for Theoretical Physics at Kyoto University. Dis- 
cussions during the YITP workshop YITP-W-04-14 on "Quantum chaos: Present 
status of theory and experiment" were useful in completing this work. This work 
is supported by a Grant-in- Aid for Scientific Research from the Ministry of Educa- 
tion, Culture, Sports, Science and Technology, Japan, for Young Scientists (B) (No. 
15740244) to H. M. 

References 

1) H. Makino and S. Tasaki, Phys. Rev. E, 67(2003),066205. 

2) M.V. Berry and M. Robnik, J. Phys. A17(1984),2413. 

3) A. Pandey, Ann. Phys. (N.Y.) 119(1979), 170. 

4) O. Bohigas, Random matrices and chaotic dynamics (LesHouches 1989 Session LIT, chaos 
and quantum physics, North- Holland) . 

5) O. Bohigas, M.J. Giannoni, and C. Schmit, Phys. Rev. Lett. 52(1984),!. 

6) M.V. Berry and M. Tabor, Proc. R. Soc. Lond. A 356(1977), 375. 

7) M.L. Mehta,_Random Matrices{2nd ed., San Diego, CA: Academic Press, 1991) 

8) F. J. Dyson and M. L. Mehta, J. Math. Phys.,4(1963),701. 

9) O. Bohigas and M.-J. Giannoni, Lecture Notes m Physics vol. 209(Springer, Berlin 1984). 

10) Ya. G. Sinai, Physica A163(1990),197. 

11) S. A. Molchanov, Commun. Math. Phys., 78,429(1981): S. A. Molchanov, Math. USSR 
Izvestija,12(1978),69. 

12) N. Minami, Prog. Theor. Phys. Suppl.ll6(1994), 359. 

13) P. M. Bleher.J. Stat. Phys.61(1990),869: P. M. Bleher,J. Stat. Phys.63(1991),261. 

14) R. D. Connors and J. P. Keating, J. Phys. A30(1997),1817. 

15) J. Marklof, Comm. Math. Phys. 199(1998), 169. 

16) J. Marklof, Progr. Math. 202, Birkhauser, Basel(2001),421; DUKE Math. J., 115(2002), 
409; Ann. of Math.l58(2003),419. 

17) A. Eskin, G.A. Margulis and S. Mozes, "Quadratic forms of signature (2,2) and eigenvalue 
spacings on rectangular 2-tori" , to appear in Ann. Math. 



12 



H. Makino and S. Tasaki 



18) 



19 

20^ 
21 

22 

23 
24 
25 
26 

27' 
28 
29 
30 
31 
32 
33 
34^ 



P.Sarnak, "Values at integers of binary quadratic forms", Harmonic Analysis and Number 
Theory (Montreal, PQ, 1996), 1810203, CMS Conf. Proc, 21, Amer. Math. Soc, Provi- 
dence, RI. 1997. 

A. I. Shnirelman, Usp. Mat. Nauk.30,(1975),265. 

G. Casati, B. V. Chirikov and I. Guarneri, Phys. Rev. Lett. 54(1985), 1350. 

V. I. Arnold and A. Avez, Ergodic Problems of Classilcal Mechaniocs, (Ben- 
jamin, N. Y.,1978): V. I. Arnold, Russ. Math. Surveys 18:6 (1963),85. 
W. Feller, An introduction to probability theory and its applications (2nd ed., John Wiley 
& Sons, Inc., New York, 1957) 

B. Li and M. Robnik, J. Phys. A27(1994),5509. 
M. V. Berry,J. Phys. A10(1977),2083. 

M. V. Berry,Phil. Trans. R. Soc. A287(1977),237. 

M. Robnik, Nonlinear Phenomena in Complex Systems l,nl(1998),l. 

A. I. Shnirelman, Usp. Math. Nauk 29(1974), 181. 
S. Zclditch, Duke. Math. J 55(1987),919. 

T. Prosen, Physica D 91(1996), 244. 

H. Makino, T. Harayama and Y. Aizawa, Prog. Theor. Phys. Suppl. 139(2000), 477. 
T.H. Seligman and J.J.M. Verbaarschot, J. Phys. A 18(1985), 2227. 

Ya. G. Sinai, Kurs teorii veroyatnostej published by MGU, Theorem 13.1 

B. V. Chirikov and D. L. Shepelyansky, Phys. Rev. Lett. 74(1995), 518. 
K. M. Prahm and D. L. Shepelyansky, Phys. Rev. Lett. 78(1997), 1440. 




Fig. 1. Numerical results for the level number variance E'^{L) with (a) Q = (1 y5)/2 and (b) 
a = 1. Wo used energy levels Cm.i G [1000000 x 10^', 1000005 x 10^]. The total numbers of levels 
are (a) 4999087 and (b) 4999224. The solid line represents the number variance in the case of 
the Poisson statistics, -S'^(i) = L. 



Long-Range Spectral Statistics of Classically Integrable Systems 



13 



c(L) 



4 
2 

-2 
-4 
-6 
-8 
-10 



Poisson statistics 



(a) 



' 00 000 0»0000000 00 oooo^ooo^oo 



(b) 



+ 

V +++++ ++ + . + + 



10 15 20 25 



Fig. 2. Numerical results for c(L) with (a) a = (1 + \/5)/2 and (b) a = 1. The solid line represents 
the number variance in the case of Poisson statistics, c(L) = 0. 




Fig. 3. Schematic picture of identification method for eigenenergy levels in the (n, m^/a) plane. 
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Fig. 4. Numerical results for c{L) with (a) q = (1 + \fZ)/2 and (b) q = 1. We used energy levels 
En,™ e [1000000 X 10^ 1000005 X lO'']. The total numbers of levels are (a) 4999087 and (b) 
4999224, and the number of spectral components is A'^ = 1000. 



